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Abstract. We have performed spatially-resolved spectroscopy of the double-nucleated Ultra-Luminous Infrared 
Galaxy IRAS 19254-7245, "the Superantennae" , along the line connecting the two nuclei. These data are analysed 
with a spectral synthesis code, to derive the star formation and extinction properties of the galaxy. The star 
formation history (SFH) of the two nuclei is similarly characterized by two different main episodes: a recent burst, 
responsible of the observed emission lines, and an older one, occurred roughly 1 Gyr ago. We tentatively associate 
this bimodal SFH with a double encounter in the dynamical history of the merger. We have complemented our 
study with a detailed analysis of the broad band spectral energy distribution of the Superantennae, built from 
published photometry, providing the separate optical-to-mm SEDs of the two nuclei. Our analysis shows that: a) 
the southern nucleus is responsible for about 80% of the total infrared luminosity of the system, b) the L-band 
luminosity in the southern nucleus is dominated by the emission from an obscured AGN, providing about 40 to 
50% of the bolometric flux between 8 and 1000 /xm; c) the northern nucleus does not show evidence for AGN 
emission and appears to be in a post-starburst phase. As for the relative strengths of the AGN and starburst 
components, we find that, while they are comparable at FIR and sub-mm wavelengths, in the radio the Sy2 
emission dominates by an order of magnitude the starburst. 

Key words, galaxies: individual(IRAS 19254-7245) - galaxies: individual(Superantennae) - interactions - starburst 
- seyfert - spectral-synthesis - SED 



^ 1. INTRODUCTION of relatively unfrequent but very luminous galaxies con- 

stitute a new important cosmological component. 



Starburst galaxies are characterized by a substantially en- . , ,., 

□ , , , • , r ,• -,i , , ,i • , With luminosities m the range 10 -10 J L n and space 

hanced ongoing star formation with respect to their past . . . ., , ° , . ,il„. 

m1 i ■ ,• n i , j i densities similar to those oi quasars (Sorter et al. 1986), 

average. 1 hough initially selected as having intense opti- , TTTT „„ , - , , \ . , 

n . . ,. ., , , , , , ,, , . the ULlRGs are the most powerful objects m the nearby 

cal emission lines, it has later become clear that the mam . it ■ — " ^ r . — r— — i . , , 

, , , r . , . j • , i Universe Sanders & Mirabel. 1996 and are considered to 

phases of star formation are highly attenuated m the op- , , , , ~r , " , ,. , , . . 

. , , , . . . „ j /,i TTi ao j- j be the local counterparts of the newly discovered high-z 

tical and luminous m the mtrarcd (the IRAS-discovered _ . . . TT , ,. , ? , 

T . , TT1 , T . T „ , „ , . Tm „ 1R luminous galaxies. Understanding the nature oi such 

Luminous and Ultra-Luminous Infrared Galaxies - LIRGs ... . . . . . , 

jttttti/^ \ rpi t -c jo ni. /to/-i\ i. ultra-luminous sources is then an important task tor o bser- 

and U LIRGs). I he Infrared Space Observatory (ISO), to- , , , , , , 

.,, , , , r ..... , „ niTrl . vational cosmology, as they may correspond to the phases 

gether with ground- based mm facilities such as SGUB A on . J , . . 

,, m ., m , , , , . r . , .,, of formation, or the assembly, ot spheroidal galaxies, 

the JGMI, showed the existence oi distant galaxies with J1 ° 

enhanced IR emission (e.g. Elbaz et al. 1999, Barger et Local infrared galaxies show evidence, from the opti- 

al. 2000, Small et al. 2000, Franceschini et al. 2001), char- cal emission lines, that their high luminosities are due to 

acterized by a strong cosmological evolution, and partly dust reprocessing of the light emitted by young massive 

responsible for the cosmic IR background (Puget et al. stars - However, the small sizes of the nuclear star forming 

1996, Hauser at al. 1998). This high-redshift population regions of compact ULIRGs are also consistent with a sig- 
nificant fraction of their IR emission being provided by an 
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Fig. 1. Position of the central slitlct on the B image taken 
at the 3.6m ESO telescope. The spatial spectral analysis 
has been obtained by cutting the slit region in 17 slices, as 
numbered in the figure. The two circular apertures centerd 
on each nuclei and wit h a diameter o f 8 arcs ecs correspond 
to those adopted by iMirabel et al See text for 

more details. 

taken as an indication for an obscured AGN contribution. 
A concomitant starburst- and AGN-activity is sometimes 
clearly revealed by broad optical and high excitation lines. 

Whether this double-nature phenomenon is occasional 
or starburst activity is systematically accompanied by an 
AGN, is still a matter of debate. The existence of a massive 
dark objects in the nuclei of galaxy spheroids (Magorrian 
et al. 1998, Merritt & Ferrarese 2001) may indeed suggest 
that the AGN activity should have been a widespread phe- 
nomenon, at least in the past. 

With the aim of further investigating the nature of 
luminous IR galaxies, we performed a spatially-resolved 
spectro-polarimetr ic study of a small sa mple of ULIRGs 
selected from the iGenzel et al. I l)l998h mid-IR sample 
of 15 ULIRGs, based on a complete IRAS sample with 
■S^opm > 5.4 Jy (Pernechclc et al., 2002). The far-infrared 
selection ensures minimal bias by dust absorption. In this 
way we aim at identifying signatures in the SEDs of these 
sources which can be used to estimate the AGN contribu- 
tion; the mechanisms quenching out the central AGN after 
the dynamical destabilization of the gas; the relationship 
between the star formation process and the dynamical in- 
teraction. 

In this paper we analyse the full-light spectrum of an 
ULIRG prototype whose nature (staburst/AGN) is still 



intriguing: the Superan tennae (IRAS 19254 -7245). At a 
redshift z = 0.061709 (|Strauss et all Il992|) and with a 
total infrared luminosity of about 1.4xl0 12 Lq (we use 
Hq = 65 [Km s _1 Mpc -1 ]), the Superantennae repre- 
sents one of the most powerful galaxy collisions in the 
nearby universe, exhibiting two giant tails extending 350 
kpc away from two nuclei, separate d by ~10 arcse c . This 
gala xy has been already studied bv IMirabel et al ] <ll99lf) 
and iGolina et al who report evidence of a type-2 

AGN with a biconical galactic wind up to 800 km/s. 

The paper is structured as follows. In Section|2]we de- 
scribe the observations and data reduction. Section [3] is 
devoted to the analysis of the most relevant characteris- 
tics of the spatially-resolved spectra and of some physical 
properties, such as extinction and star formation rate. In 
Section|3]we build a consistent broad band spectral energy 
distribution of the two nuclei. These allow us to evaluate 
in Section El the contribution of the star formation pro- 
cess and the AGN to the bolomctric infrared emission. 
The optical spectrum of the two nuclei, corrected for the 
contribution of the AGN, is thus analyzed by means of a 
population synthesis tool described in Section |5j the re- 
sults are presented in Section [7| Section [5] is devoted to 
the discussion of the star formation history and the AGN 
non-thermal component of Superantennae. Section [9] con- 
tains our conclusions. We adopt H n = 65 [km s _1 Mpc -1 ]. 

2. OBSERVATIONS AND DATA REDUCTION 

The spectropolarimetric observations of the 
Superantennae were done with EFOSC2 at the 3.6m ESO 
telescope in La Silla, during the night of November the 
1st, 2000. The seeing was about 1 arcsec and some veils 
were present at low sky altitude. For this reason partic- 
ular care has been taken in observing the target at its 
maximum elevation and chosing the standard star close to 
the object, in order to perform absolute flux calibration. 
We notice in advance that, once properly corrected for 
aperture losses, the calibrated spectra consistently fit the 
available literature optical and near-IR photometry (see 
sect. El- 

Polarimetric observations were obtained adopting a ro- 
tating A/2 retarder plate and a Wollastone prism, inserted 
before the grism. We choose the 236 lines mm -1 grism, 
with a dispersion of 2.77 A pixel -1 , covering the 3685- 
9315 Awavelength range and a 1.2x20 arcsec slit, corre- 
sponding to a spectral resolution of 2 1.2 A. The slit was ori- 
ented with a P. A. of 14°, such that both Superantennae's 
nuclei were included. 

The data consist in a set of four couples of frames, 
obtained with retarder-plate angles of 0, 22.5, 45 and 
67.5 degrees, for a total integration time of 6000 seconds. 
Johnson B and R band images were also obtained, with 
exposure times of 900 and 480 seconds, respectively. 

The data reduction has been carried out with the use 
of iraf's standard tasks and idl procedures developed by 
ourselves. 
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Fig. 2. The three-dimensional structure of the spectrum obtained with the pipeline described in Section Different 
emission lines are already clearly visible, Ha standing over all the others. The spatial axis corresponds to the one 
shown in Figure ^ 



Polarimetry is reported by IPernechele et alJ l|2002() . 
while here we analyze the full-light spatially-resolved spec- 
trum. 

The Wollastone prism splits the light into two perpen- 
dicularly polarized beams, which are projected onto the 
CCD. To avoid overlapping, the slit is transparent only in 
a series of 5 slots, separated by 5 opaque masks. The prism 
causes each frame to appear as composed of 5 couples of 
spectra (hereafter "strips"), belonging to the 5 different 
slots. 

The frames were bias and flat-field corrected in the 
usual way. Then the two strips including the object's 
traces and two sky-strips were cropped. Each strip was 
flux calibrated with the spectr opho t ometr ic standard star 
in the centre of NGC 7293 l)Okel 11990^) and then sky- 
subtracted. Unfortunately, the interposition of the prism 
between the retarder plate and the grism induces strong 
distortions of the frames, so we needed to provide a dis- 
tortion map by means of a polynomial fit using idl. Each 
spectrum was corrected for gal actic reddening, assu ming 
Rv = 3.1 and E(B-V)=0.086 l|Schletrel et al.Ul99£fr and 
finally the various strips, relative to 4 different polariza- 



tion angles, were combined together to produce the spa- 
tially resolved total- light spectrum of Figure [3 

The spectra of the nuclei were simply extracted 
through apertures of 1.2x3.6 arcsec 2 (for the southern) 
and 1.2x2.4 arcsec 2 (for the northern nucleus), see Figure 
03 At a redshift z = 0.06171, 1 arcsec corresponds to ~ 1.4 
kpc. 

3. ANALYSIS OF THE SPATIAL DEPENDENCE 
OF THE MOST RELEVANT SPECTRAL 
FEATURES 

The high quality total-light spectrum allowed us to 
perform an analysis of the main spectral features of 
Superantennae along the whole slit, hence to probe the 
spatial dependence of some physical properties along the 
line connecting the two nuclei. 

To take into account the effects of the poor seeing, we 
binned our spectrum in steps of 3 pixels along the spa- 
tial axis (corresponding to 1.2 arcsec) and thus obtained 
seventeen independent spectra corresponding to different 
regions of the galaxy. In Figure Q the 1.2" slit has been 
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Fig. 4. The seventeen different slices of Figure Q are plotted in sequence, to show the trend of the continuum and 
other spectral features, as a function of position along the slit. 



superposed on the B band image taken with EFOSC2: the 
different spectra obtained with the binning procedure cor- 
respond to the regions numbered from 1 to 1 7 in the figure. 
This a lso indicates the 8" apertures used bv lMirabel et alJ 
in their photometric study of the Superantennae 
(which we will use in the following to estimate the contri- 
bution of the two nuclei to the total infrared flux) . 



3.1. Emission lines 

The 17 spectra corresponding to the different regions along 
the slit (FigUJ are shown in Figure^] Several optical emis- 
sion lines typical of starburst galaxies can be recognized 
(see also FigEJ: [On] (3727 A), H/3(4861 A), the sum of 
the two [Om](4959 and 5007 A), [Ol](6300 A), the blend 
of Ha (6563 A) and the two [Nn]'s (6548 and 6584 A), 
and the blend of [Sn] (6717 and 6731 A). 

The corresponding line fluxes are reported in Table 1. 
Figure compares the profiles of the various line fluxes 
(solid line) along the spatial direction with the intensity 



of the continuum between 5500 and 5600 A (dotted line) 
normalized to the line intensity peak value. It turns out 
that, with the exception of the [On] (see Fritz et al., 2002, 
for a discussion), all other emission lines are much more 
peaked than the continuum, indicating that the ongoing 
star formation or the AGN activity is strongly concen- 
trated within the nuclei. 

The southern nucleus shows both narrow permitted 
and forbidden lines; we don't have enough spectral reso- 
lution to separat e Ha (6548 A) from t he two [Nn]'s (6563 
and 6584 A), but lMirabel et alJ l|l99l|) reports that [Nil] is 
stronger and sharper than Ha, thus revealing the Seyfert 
2 nature of this emission. 

For the same reason we cannot discriminate between 
[Om]A4363 and H 7 or [Sn]A6717 and [Sn]A6731 and, 
consequently, we cannot derive an estimate of the elec- 
tron density. By adopting the va lue of iV e = 670 cm~ 3 
provided by l|Colina et all Il99l|) . we use the fluxes of 
[Nn]A5755 and [Nil] AA6548+6583 to estimate the elec- 
tron density in the southern nucleus (spectra #4, 5 and 
6). [Nil] A5755 was deblended from [Fevn] A5721, while the 
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Fig. 5. Spatial profile of the main emission lines. The solid histograms represent the profiles of the various lines, as 
indicated in the panels. To study the spatial profile of the lines, we plot the continuum normalized to the peak value 
of line emission in the sounthern nucleus (dotted histogram): going from the brightest nucleus to the outer regions, all 
the lines but [On] fade more rapidly than the continuum. 



ratio between Ha and [N il] AA6548+6583 has been inferred 
from lColina et af] l)l99lh data as described next. The es- 
timated electron temperature in the southern nucleus is 
T e ~ 14600 K, consistent with the Colina et al. (1991) 
range of 10 4 - 2.5 10 4 K. 

The northern nucleus is characterized by a much lower 
star formation activity, as it is evident from the lower in- 
tensity of the lines or even the prevalence of photospheric 
absorption for H/3 and H7. A detailed analysis reveals that 
the [Oni]A4363 + H7 blend is in emission in the three cen- 
tral spectra of the southern nucleus, possibly because of 
a strong forbidden oxygen emission in the Seyfert compo- 
nent, but then moving outwards and to the northern nu- 
cleus it switches to absorption, dominated by photospheric 
H7. Also the blue edge of the continuum steepens signifi- 
cantly towards the northern component, showing that this 
nucleus is dominated by an older stellar population. 

3.1.1. Extinction from the line ratios 

The comparison of the observed Ha/H/3 ratio with the 
theoretical value of 2.85 predicted for a Case-B recombi- 



nation l|Osterbrocklll989|) . allows in principle to derive the 
average extinction, at least for the southern nucleus where 
H/3 is detected in emission. To this end, we first need to 
quantify the contribution of the two [Nil] lines at 6548 
and 6584 A to the Ha+[Nil] blend, because the Seyfert 
2 nature of this nucleus prevents from adopting standard 
starburst values for the line ratio. F rom the high resolu- 
tion spectrum of lColina et alJ l)l99l|) , we know that in the 
southern nucleus 5(Ha) = |5([Nn](6584)), hence assum- 
ing for the ratio between the two [Nil] the the standard 
value (|), we get: 

S(Ea) = — S(blend). (1) 
25 

Simil arly, in the northern nucleus we obtain l|Colina et all 
Il99lh S(Ha) =s |5(blend). 

Since Colina et al. (1991) spectra of the two nuclei were 
obtained through a 3"x6" aperture, in the following we 
adopt eq.Q]on spectra #2 to #8, centered on the Southern 
nucleus and the more standard value of h in the Northern 
nucleus and outer regions of the Superantennae. 

To compare the empirical extinction estimate 
with the model prediction (see Section |8J we need 
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SfB~) 










SfHa+fNnl") 




E(B-V) 


1 


6.68 


0.28 




2 01 +2 84 


0.36 




15.3 


4.28 




2 


5.88 


0.34 


1.28 


2 47+4 69 


0.43 


1.53 


29.3 


8.03 


0.429 


3 


11.1 


0.41 


2.28 


3 32+8 16 


0.55 


4.97 


72.9 


18.0 


0.730 


4 


14.4 


0.57 


7.65 


16 1+36 7 


0.88 


25.8 


295. 


57.9 


1.013 


5 


19.3 


0.85 


20.5 


42 3+99 3 


1.47 


74.5 


826. 


125. 


1.090 


6 


16.8 


0.75 


16.5 


28 4+72 7 


1.25 


46.8 


570. 


87.8 


0.974 


7 


11.6 


0.43 


4.10 


7 1 n+ 1 8 5 


0.68 


7.13 


130. 


22.4 


0.827 


g 


77.0 


0.31 


1.69 


2 70+5 60 


0.49 


3.20 


36.7 


7.37 


0.450 


g 


6.46 


0.28 




n 814+2 91 


0.44 


1.31 


17.3 


4.48 




10 


5.21 


0.28 




1 03+2 78 


0.42 


1.06 


14.0 


3.73 




11 


5.04 


0.29 






0.41 




13.3 


3.53 




12 


4.91 


0.34 






0.49 


1.55 


14.5 


4.09 




13 


5.15 


0.53 


-1.50 


3.70+4.51 


0.82 


1.75 


26.3 


7.54 




14 


5.52 


0.66 


-4.20 


2.17+2.87 


1.09 


1.31 


33.5 


8.26 


- 


15 


5.58 


0.56 


-2.80 


n 885+1 36 


0.73 




16.9 


6.12 




16 


5.75 


0.42 






0.47 




9.44 


3.13 




17 


5.51 


0.30 






0.39 




9.52 


2.84 




South. Nucl. 
North. Nucl. 


48.8 
9.26 


2.16 
1.19 


43.9 
-6.60 


84.8+208 
3.73+5.03 


3.59 
1.91 


156. 
4.03 


1680. 
57.0 


240. 
13.2 


1.042 



Table 1. Flux intensities of the main emission lines detected, measured on the extracted spectra (no slit losses 
correction applied yet). The first column indicates the label of the spectrum as specified in Figure ^ I n the next 
eight columns we report the observed fluxes of [On] (3727 A), Rf3 (4861 A), [Om] (4959 + 5007 A), [Oi] (6300 A), 
Hev+[Nn] (6563, 6548 and 6584 A) and [Sn] (6717 A), respectively. Fluxes are in 10~ 16 [erg cm" 2 s -1 ]. The last 
column shows the E(B-V) derived from the Balmer decrement, as described in the text. The last two rows report the 
fluxes for the two nuclei: the southern consists of spectra #4, 5 and 6, the northern of # 13 and 14. The nuclear fluxes 
are larger than the sum of fluxes in individual spectra, because they were measured on total nuclei extractions and 
not computed as sums. 



Southern Nucleus h«+[nii] 




4000 5000 6000 7000 8000 

A(A) 



Fig. 3. The spectra of the two nuclei of IRAS 19254-7245 
as extracted from an aperture of 1.2" x 3.6" for the south- 
ern nucleus, and 1.2" x 2.4" for the northern, respectively. 
A lot of emission lines typical of ULIRGs are recognizable. 



to further correct the observed H/3 flux for the 
photospheric stellar absorption. On this subject 
lOonzalez Delgado. Leitherer. fc Heckmanl l|l999ft found 



that the contribution to the H/3 equivalent width from 
intermediate-age stars in starburst galaxies is 2 — 5 A. 
Since they do not account for an old disc population, 
which would decrease these values, we have adopted a 
mean absorbing equivalent width of 2 A. The color excess 
E(B-V), derived from the Hev/H/3 ratio is reported in the 
last column of Table 1. We plot in Figure the spatial 
dependence of the E(B-V) (upper panel) and of Ha 
(middle panel) and H/3 line intensities. 

The values of E(B-V)~1.0 mag in the central regions 
obtained in this way ar e smaller than the va lue of E(B-V) 
~1.4 mag reported by IColina et alJ lll99lh . but are con- 
sistent with those obtained by our more detailed spectral 
analysis discussed in Sect. w hich conside r s the pres- 
ence of an old population. Also IVanzi et al.l l|2002h find 
Ay = 3.15, consistent with our result. This technique can- 
not be applied in the outer regions of the system and in 
the northern nucleus, where H/? is either not detected or 
in absorption. In those regions, spectral synthesis models 
(see fig. E3> yield a mean value of E(B-V)~ 0.4 — 0.5, 
consistent with the typical average colour excess ~0.3 of 
spiral galaxies. 

Hopkins et al (2001) and Sullivan et al. (2001) notice 
the existence of a correlation between the extinction, de- 
rived from the Balmer decrement, and the intensity of the 
Hcv emission line, in a sample of star forming galaxies. 
The same trend is clear in our spatially resolved analysis 
of Superantennae, when we plot the ratios of [On] (3727A) 
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-6-3 3 
arcsec 



Fig. 6. Upper panel: spatial dependence of the extinction, 
E(B-V), as computed from the Balmer decrement along 
the slit. Middle and lower panels: Ha and H/3 fluxes. The 
H/3 flux scale is linear instead of logarithmic because the 
line is seen in absorption in the Northern nucleus. 



and H/3 to Ha, as a function of the Ha flux (Figured). 
The observed ratio decreases as the Ha intensity increases, 
confirming a strong correlation between extinction and the 
star formation activity, as expected if the latter is a func- 
tion of the local density of gas and dust, and if the extinc- 
tion is mainly due to a foreground screen. The correlation 
involving [On] (3727A), see top panel in Fig. [7| seems to 
indicate that [On] and Ha emissions are largely indepen- 
dent, in fact the [On] strength only slightly depends on 
position along the slit. We defer to Fritz et al. (2003) for 
a more detailed discussion of this effect. 

3.1.2. The Star Formation Rate 

In star forming galaxies the extinction-corrected Ha line 
intensity is a powerful tool to estimate the rate of ongoing 
star formation. Poggianti, Bressan & Franceschini (2001), 
however, have warned that for particularly extinguished 
objects the star formation rate determined in such way 
may be significantly underestimated, by values up to a 
factor 3 compared with those derived from the far infrared 
luminosity. Furthermore, as it will be discussed later, the 
Superantennae shows strong evidence of a significant AGN 
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Fig. 7. [Oll](3727A)/Ha and H/3/Ha flux ratios vs. Ha, 
before extinction correction. Each data point is labelled 
according to Figured (the southern nucleus refers to spec- 
tra #4, 5, 6; the northern to # 13 and 14). There is a clear 
dependence of the ratios to the intensity of Ha, suggesting 
a tight relationship between extinction and Ha luminos- 
ity, indicative that the most active (luminous) regions are 
also the more extinguished. See Fritz et al. (2003) for a 
detailed discussion. 

component, probably contributing also to optical emission 
lines. 

Wit h these caveats i n mind, we have applied the rela- 
tion of lKennicuttl l|l998j) to the Ha emission corrected for 
the [Nil] contribution to the blend, for extinction, and for 
the AGN contribution, to derive a preliminary estimate of 
the corresponding SFR, leaving a more accurate analysis 
to the following Sections. 

Kennicutt assumes a Salpeter IMF between 0.1 and 
100 M Q , solar abundance and a constant star formation 
rate during 10 8 yr: 

SFRiM^yr- 1 ) = 7.9 • lO" 42 ^ (ergs" 1 ). (2) 

The results are summarized in the first column of Tabled 
(see Section |8J, with values ranging between ~ 47 M & /yr 
in the southern nucleus and ~ 0.4 M^/yr in the northern. 

4. THE BROAD BAND SPECTRAL ENERGY 
DISTRIBUTION 

Further information on the source comes from the analysis 
of the broad band energy distribution. We have collected 
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all the photometric data available on IRAS 19254-7245 
with the purpose of obtaining a reliable spectral energy 
distribution for its two components. 

IMirabel et all lll99lft report magnitudes in the opti- 
cal and near-IR, from the B up to the L and N bands, 
within 8" apertures centered on the two nuclei, as well as 
the total corrected IRAS fluxes at 12, 25, 60 and 100 fim 
with uncert a inties of 11%, 5%, 4% and 8%, respectively. 
iKlaas et alJ l|200ll) publish ISOPHOT measurements for 
the whole object at 10, 12, 15, 25, 60, 90, 120, 150, 180 
and 200 /im and SEST data at 1300 /mi. They quote un- 
certainties of about 5% at 100-150 fim and up to 10%-20% 
for the other ISOPHOT passb ands and for SEST. 

ICharmandaris et all ()2002f) report ISOCAM observa- 
tions in various mid-IR bands of the interacting system 
with a spatial resolution of few arcsec, hence resolving 
the two nuclei. These results are consistent with the pho- 
tometric observations in the N-band by Mirabel et al. for 
the southern nucleus, but are systematically fainter for the 
northern one. Then we split the IRAS and the ISOPHOT 
long wavelength fluxes between the two nuclei based on 
the Charmandaris et al. photometry: for A < 12 /im we 
assumed a ratio between the IRAS and ISOPHOT fluxes 
belo nging to the two nuclei e qual to the ratio between 
the ICharmandaris et all l)2002(l ISOCAM 11.4 /xm fluxes; 
while longward 12 /im we ad opted the ratio between the 
fluxes from the two nuclei bv ICharmandaris et alJ ((2002) 
in the ISOCAM 15 /im filter (see tableland Fig. EH). 

Since the northern component is very faint compared 
to the southern one, the uncertainties in method strongly 
affect the fluxes of the northern nucleus. Therefore, as 
concerns the northern nucleus, we consider fluxes obtained 
in this way reliable only shortward of 25 /Ltm. 

Almost 90% of the infrared flux is emitted within the 
two 8" circular regions centered on the nuclei to which the 
optical photometry is referred, the southern nucleus being 
responsible for about ~ 95% of this flux. 

4.1. Slit-loss corrections 

To match the broad band SED obtained above with the 
optical spectra measured within our slit, we have mod- 
elled the surface brightness profile of Ha with an axially- 
symmetric bidimensional gaussian. The free parameters, 
namely the peak value and the variance er, were estimated 
by fitting this profile along the slit with a linear gaussian 
as shown in Figure |H1 The values derived are a sou th = 2.6 
pixels and a nor th = 2.1 pixels, corresponding to 1.04" and 
0.84", for the southern and northern nuclei respectively. 

The volume under the 3D gaussian is simply given by 
the integral of the profile summed to the zero-level flux 
given by the continuum between the nuclei. To evaluate 
what fraction of the emitted flux falls outside the slit, we 
integrated "by slices" within the 1.2" x 8" apertures cen- 
tered on the nuclei. As a result, about the 58% of the 
light coming from the 8" aperture centered on the south- 
ern nucleus, and about 75% of that coming from the same 
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5.205-10~ l;i 


1.26-10~ 14 
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Table 2. Observed fluxes and magnitudes of the 
Superantenna e. Op tical and IRAS data, from 
IMirabel et alJ l|l99lt) . are in the first thirteen rows, 
the next six rows repor t the IS OCAM LW bands fluxes 
bv ICharmandaris et alJ l|2002() and finally IKlaas et alJ 
(200l|) measurements are in the last eleven rows; the 1300 
/xm datum was obtained with SEST, while the others 
with ISOPHOT. Infrared observations marked with * 
have been splitted between the two nuclei, by means of a 
procedure described in the text. Integ rated IR fluxes were 
obtain ed from IRAS data through ISanders fc Mirabel 
l(l996l) calibration and are expressed in W-m -2 . 



region around the northern nucleus, were lost by our op- 
tical spectroscopic observations. 

The next step was to consider the amount of light en- 
closed in the observed nuclear regions, namely within spec- 
tra numbered 4, 5 and 6 for the southern component, and 
13, 14 for the northern one. Considering that the observed 
Ha emission is likely to appear broadened because of the 
higher extinction in the very inner region of the nucleus, 
we find that about 50% of the flux observed from a circu- 
lar aperture of 8" on the southern object is emitted within 
this 1.2" x 3.6" aperture corresponding to spectra 4, 5, and 
6. Assuming that the far-IR emission mimics the spatial 
distribution the Ha line, the ratio S1R/S5550 between the 
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Fig. 8. Spatial profile along the slit, at the wavelength 
of Ha, used to evaluate the amount of energy encircled 
within the 8" aperture centered on each nucleus (delim- 
ited with shaded borders), which falls outside the slit (1.2" 
shaded regions). After performing a gaussian fit to the 
profile, slit losses were estimated assuming that the bidi- 
mensional light profile is an axial-simmetric gaussian. 



bolometric infrared and the V band fluxes turn out to be 
~ 90 for the southern nucleus 1 . 

Similarly we estimate a slit-loss of about 90% from the 
northern nucleus (1.2" x 2.4" aperture for spectra 13 and 
14) and, correspondingly, a ratio S1R/S5550 ~9. 

5. THE SPECTRAL SYNTHESIS MODEL 

The optical spectrum carries a wealth of information on 
the recent star formation history and dust extinction. In 
this Section we describe our detailed analysis of the optical 
spectrum for each separate nuclear component. 

As a best approach to describe the irregular star- 
formation histories of starburst galaxies, we have 
adopted a completely free-form galaxy spectral synthe- 
sis model derived as a modofication of that discussed in 
IPoggianti et alJ l)200l|) . The optical observed spectrum is 
modelled as a combination of up to 10 simple stellar pop- 
ulations (SSP) of solar metallicity and different age. Each 
SSP is meant to represent a formation episode of average 
constant star formation rate (SFR) over a suitable period 
At (see Table 0). 

As the main information provided by the spectrum are 
the shape of the continuum, (nebular) emission lines and 
(photosphcric) absorption lines, we have chosen 3 classes 
of SSPs spectra to reproduce these features: young, in- 

1 The flux at 5550 A is computed as the mean value within 
about 50 A centered at A = 5550 



termediate and old. The ages of the 10 populations have 
been chosen by considering the evolutionary timescales 
of the stars associated with the observational constraints: 
the youngest generations (10 6 , 3 • 10 6 , 8- 10 6 , 10 7 yr) are re- 
sponsible for the ionizing photons producing the emission 
lines an for the UV continuum; intermediate populations 
(5-10 7 , 10 8 , 3-10 8 , 5-10 8 , 10 9 yr) have the strongest Balmer 
absorption lines; while older generations of stars provide a 
significant contribution to the spectral continuum, hence 
affecting also the equivalent widths of the lines. The lat- 
ters have been modelled with a constant star formation 
rate (SFR) between 2 and 12 Gyr before the observation. 

Each single population is assumed to be extinguished 
by dust in a uniform screen according to the standard ex- 
tinction law of the diffuse medium in our Galaxy [Ry = 
A V /E(B - V) = 3.1, Cardelli et al. 1989]. While a more 
complex picture of the extinction cannot be excluded, 
Poggianti et al. have already shown that, in the case of 
obscured starbursts, the characteristics of the emerging 
spectrum require a significant amount of foreground dust 
(screen model). For a uniform mixture of dust and stars, 
an increase of the obscuration does not yield a correspond- 
ing increase in the reddening of the spectrum, and E(B-V) 
saturates to a value of ~ 0.18, too low to account for the 
observed emission line ratios. In our case the extinction 
value E(B-V) is allowed to vary from one stellar popu- 
lation to another and the summed extinguished spectral 
energy distributions of all the single generations make up 
the total integrated spectrum. 

The integrated spectra of the SSPs have been com- 
puted with a Salpeter initial mass function (IMF) be- 
tween 0.15 and 120 Mq. We have adopted the Pickles 
(1998) spectral library, after extending its atmospheres 
with Kurucz (1993) models outside its original range of 
wavelengths, from 1150 A to 25000 A. The composite 
(stars+gas) spectra have been obtained through the ion- 
ization code CLOUDY (Ferland, 1990). 

The best-fit model was obtained by minimizing the 
differences between observed and model selected features: 
the equivalent widths of five lines ([OII]A3727,H5, H/3 and 
Ha) and the relative intensities of the continuum flux in a 
number of almost featureless spectral windows which are 
chosen according to the spectra charachteristics. 

As already anticipated our spectr al resolution does no t 
allow to resolve the Ha- [Nil] blend. IColina et all l)l99l|) . 
on high resolution spectra, have found that Ha is about 
30% fainter than [Nil] 6584. Because our spectral synthe- 
sis analysis does not model the [Nil] 6548 line, we apply 
a correction factor to account for this contribution to the 
blend. 

A critical constraint that we impose to the model is 
to reproduce the observed infrared to optical flux ratio 
S1B./S5550. Our predicted infrared flux is computed as 
the difference between the total non-absorbed and the ab- 
sorbed model spectra. 

The minimization procedure appl ies the algorithm o f 
Adaptive Simulated Annealing (ASA. Ilngberl Il989t 120001 
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Fig . 9. Spectral energy distribution of the southern nucleus. Data are from lMirabel et all 1 199lh . [Laurent et alJ |2000) 
and lKlaas et alJ l|200ll) as described in the text. The thick solid line is the fit to the global SED, obtained combining 
an optical-NIR fit (dashed line short ward 5 /xm) to our spectra (corrected to the 8" aperture) with an AGN model 
obtained with the code DUSTY (dot-dashed line) and a semiempirical starburst Arp220-like template (dashed line 
long ward 5 /im). The thin solid and dotted lines show the fit of the central and outer regions spectra. The N band 
width is slightly shifted upward with respect to the observed datapoint, in order not to be confused with the LW2's. 



to a suitable merit function (MF) based on the differences 
between modelled and observed selected spectral features: 

where Mj and Oi are respectively model and observed 
fluxes (or equivalent width), and Ei and Wi are the cor- 
responding errors and weights. The latter are used to re- 
duce the weight of some observables particularly difficult 
to model (as the Oil line flux). 

The minimization procedure checks if any of the SSPs 
older than 10 7 years contibutes by less than 1% to the 
total flux at each wavelength. This may happen, for in- 
stance, when extinction becomes too large; in that case 
that SSP is no more relevant for the fit and is not used. 
Obviously, the same does not apply to younger SSPs be- 



cause, though very extincted, their contribution to the 
FIR emission could not be neglected. 

Effectively only few populations dominate the sinthctic 
resulting spectrum: as can be seen from table [21 the num- 
ber of SSP used in the fits decreases to 4 or 5 tipically. 



6. THE AGN CONTRIBUTION TO THE BROAD 
BAND SPECTRAL ENERGY DISTRIBUTION 

The broad band SEDs of the two nuclei of the 
Superantennae, together with the best-fit optical synthetic 
spectra are shown in Figures 1^1 and ITU1 respectively. 



Berta S. et al.: Optical Spectroscopy of the Superantennae 



11 



6.1. The southern nucleus 

The resulting SED of the southern nucleus is compared in 
Figure IH1 with a model obtained from the combination of 
three different physical components. 

In the optical domain, it is the result of our spectral 
synthesis analysis. The model is obtained combining two 
different fits, belonging to the central (# 4, 5, and 6 - 
thin solid line in Figure |5J) and outer (# 2, 3, 7 and 8 
- dotted line) regions spectra, normalized in order to ac- 
count for the slit losses. The two spectral combinations 
are used to describe the nuclear and large-scale galactic 
emissions. Indeed the sum of the two (short-long dashed 
line) reproduces quite well the optical photometric data, 
as anticipated in section |2J It is worth noticing that the 
spectral model reproduces well also the photometric data 
beyond 7000 A up to the H band. 

A first inspection shows that from this fit there re- 
mains a strong residual excess in the K and L bands with 
respect to the stellar component fitting the optical data. 
Longwards of the H band, the observed fluxes seem to 
follow a power-law shape, possibly due to the Seyfert II 
emission, while the stellar spectral component turns down. 
This excess is even more evident in the L band, where 
the stellar emission is fainter and dust reprocessing of 
young stellar emission is not yet relevant. The observed 
J-L colour is about 3.9 mag, while that of the best fit galac- 
tic model is J-L~2.7mag. The latter is consistent with the 
value observed in the northern nucleus, see Table and 
Fig. which does not require an AGN contribution. 

In order to reproduce this excess we added a QSO 
template, requiring that the stellar plus AGN compo- 
nent, due to circumnuc lear dust, match the L band flux. 
IPernechele et alJ l|2002() performed a detailed polarimetric 
analysis of the optical emission of IRAS 19254-7245, re- 
vealing a scattered polarized component at 2% level only. 
For this reason we generate a totally obscured AGN model 
with the code DUSTY (Ivezic et al., 1999) assuming a 
spherical geometry and a standard ISM. As an alterna- 
tive, we also checked our fit using the totally obscured 
edge-on AGN models by Andreani et al. (1999), finding a 
full consistence with the DUSTY results. 

The free parameters for the AGN fit are the ratio R 
between the inner and outer radii of the spherical dust 
distribution and the optical depth t\ at 0.55 /an. 

At wavelengths longer than 10 /zm the contribution 
of the starburst starts to rise again. We thus added an 
Arp220-like starburst template (short-long dashed line) 
normalized in order to account for the remaining IR lu- 
minosity produced by the southern nucleus. Our template 
was obtained by m erging the GRASIL synthetic model by 
I Silva et all Il998|) of Ar p220 to the CAM-CVF spec trum 
between 5 and 20 /urn bv lCharmandaris et alJ (^99). 

We then imposed that the sum of the AGN and star- 
burst templates should reproduce the observed optical- 
infrared SED. Our best fits for the AGN emission were 
obtained for R values between 150 and 400, and an ab- 
sorption To.55 M m = 30 — 40 mag. Adopting higher values of 



To.55/im would underestimate progressively the K band and 
the L band fluxes, as r grows, because of self absorption 
in the AGN dust emission. By converse, lower extinction 
values would lead to an overestimation of the NIR fluxes. 

The parameter R, modelling the size of the circum- 
nuclear dust distribution around the AGN, has a moder- 
ate impact on our global fitting procedure: by increasing 
R we essentially add cold dust in the outer parts far away 
from the central AGN, hence slightly increase the AGN 
emission around 100 /an. However, there is limited vari- 
ance allowed to this: our best fit shows that with values 
R<150 or R>400 would completely spoil the fit to the 
observed SED. 

In this way we estimate that the contribution of the 
AGN to the bolometric IR flux (estimated from IRAS 
data) ranges between about 40% and 50%. 

The resulting model is the thick solid line in Figure El 
we note that it does not fit the observed 6—7 /im data, 
possibly characterized by strong PAH features. Actually 
DUSTY doesn't account for the PAH emission and there 
could be also some residual PAH from the starburst com- 
ponent, not reproduced by the Arp220 template. On the 
other hand the cause of this discrepancy could also be the 
spherical geometry assumption. 

This analysis, in turn, suggests that a similar fraction 
of the Ha emission could originate from the AGN. Indeed 
Cid Fernandes et al (2001), analysing the properties of a 
sample of 35 Seyfert II nuclei, find that it is possible to 
set some rough constraints on the non-thermal H/3 com- 
ponent from the properties of the [Om]/H/? ratio: they 
show how values of [Oiii]/H/3 between 4 and 13 imply an 
AGN contribution to H/3 between 30 and 100%. On this 
basis, the measured IRAS19254 South [Om]/H/3 of - 6 
leads to a non-thermal contribution of ~ 40% to the H/3 
emission, consistent with our estimate in the infrared. 

In the followings we will thus assume that about 40% 
of the line emission is due to the AGN. 

6.2. The northern nucleus 

As far as the northern nucleus is concerned, we found es- 
sentially no evidence for an AGN component. The ob- 
served data consist in the set of optical and near-IR ob- 
servations by Mirabel et al. (1991), the mid-IR ISOCAM 
fluxes by Charmandaris et al. (2001) and the IRAS and 
ISOPHOT fluxes up to 25 /im, obtained splitting litera- 
ture data for the whole (North + South) Superantennae 
through the method described in section 0] The latters 
are obviously more uncertain, but are still consistent with 
Charmandaris et al. (2001) data, which are affected by a 
~ 10 — 20% uncertainty. 

The fit is obtained by means of a pure stellar compo- 
nent, modelled as in Section|5l from the optical to the NIR, 
and a starburst component based on the prototype M82, 
at A > 5/im. The lat ter is obtained by merging a synthetic 
SED modelled with ISilva et alJ l|l998h . with the observed 
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Fig. 10. The spectral energy distribution of the northern 
nucleus. The fit (thick solid line) is obtained by adding an 
optical-NIR model short ward 5 /im, and a semiempirical 
starburst M82-like template (long ward 5 /xm). The thin 
solid and dotted lines show the fit of the central and outer 
regions spectra. No AGN component is needed in this case. 



CVF spectrum bv lSchreiber et all l|200lj) between 5 and 
18 /jni. 

The best fit is shown in Figure EH A small excess still 
present in the L band may be due to PAH 3.3 /mi emission 
not included in our spectral model. 

7. RESULTS OF THE SPECTRAL ANALYSIS 

In this Section we discuss our results for the two nuclei 
separately and then we extend our analysis to the outer 
regions of the galaxy. 

The merit function (equation |3J) has been built 
on 10 continuum bands (3550-3670 A, 3890-4020 A, 
4140-4300 A, 4400-4500 A, 4700-4820 A, 5080-5180 
A, 5360-5660 A, 5950-6110 A, 6880-6980 A and 
7750-8500 A), and 5 spectral lines, Ha+[Nll], H/3 and 
[On] 3727, H<5 and H7. The latter two lines were assigned 
a weight of 0.3 because they are too faint and, in the case 
of H7, because of contamination by [Oiii]4363, which is 
not considered in our models. 

We required that the model should reproduce the ratio 
Sir/ S5550 between the far-infrared and the optical flux, 
as defined in Sect. 0] and aperture corrected according to 
Sect. 14. ll The FIR flux was corrected for the contribution 
of the AGN, assumed to amount to 40% of the bolometric 
IR emission (see Sect.[f)J). 

The errors in the continuum flux were computed ac- 
cording to the local signal to noise ratio, while in the case 
of the lines fluxes the relative errors were assumed pro- 



portional to l/y^A, where I a is the flux emitted by the 
corresponding line. 

Table[21summarizes the parameters leading to the best 
fits: for each spectrum SFR and Ay of the populations 
which build the sinthetic model are reported. Basically 
only 4 or 5 SSPs significantly contribute to the emitted 
spectra, therefore decreasing to 8 or 10 the number of 
parameters used to fit the 16 observed features chosen. 

7.1. Models for the southern nucleus 

The spectrum of the southern nucleus is obtained by sum- 
ming the spectra with the strongest continuum and emis- 
sion lines, namely those labelled with numbers 4, 5 and 6 
(see Fig.[TJ). 

Our best-fit solution is illustrated in Fig.^] The upper 
panel shows the observed and model spectrum (both nor- 
malized at 5500 A) and their differences, the lower panel 
reports the contribution of the dominant stellar popula- 
tions. All the modelled features are within ler error from 
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Fig. 11. Best fit model for the southern nucleus. Top panel 
reports the observed (thick line) and the model (thin line) 
spectrum, the dotted line shows the residuals between 
model and observed spectra. Bottom panel depicts the 
contribution of each extinguished SSP (labeled from 1 to 
10 with increasing age). The intermediate age population 
dominates the continuum luminosity, while the 2-6 Myr 
population provides the emission lines. 

their observed values, except for the [Oil] line (3cr) and 
for the band around 4760 A. Fig. ^| illustrates the star- 
formation history and extinction distribution of the best- 
fit model. The optical spectrum is mainly contributed by 
an old population (about 12 Gyr), providing about 65% of 
the total luminous mass, and an intermediate one (about 
5 • 10 s — 10 9 yr) forming about 35% of the total lumi- 
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Fig. 12. SFR and extinction of the southern best fit 
model. Top panel: SFR (histogram and left scale) and 
colour excess (filled squares and right scale) for the differ- 
ent stellar populations. Middle panel: ratio between the 
bolometric FIR emission due to each SSP and the total V 
band flux. Notice that the sum over all populations is not 
100, but turns out to be 57.6, which is the ratio between 
the total FIR and V fluxes. This number is the conversion 
factor to obtain the differential contribution of each SSPs 
to the integrated infrared light. Bottom panel: cumula- 
tive mass formed in the galaxy: about 35% of the total 
stellar mass has been processed by the intermediate age 
star formation episode, 60% belongs to the oldest. The 
youngest stars, significantly contributing to the FIR bolo- 
metric emission and dominating the lines fluxes, provide 
almost no contribution to the total luminous mass formed 
in the galaxy. 



Fig. 13. A good fit to the optical spectrum of the south- 
ern nucleus (not shown here) is found even by imposing 
a constant long burst of star formation. However the fit- 
ting procedure "obscures" unnecessary populations, which 
only contribute to the FIR. Due to this behaviour of the 
color excess we exclude such a star formation scenario. 



value, and there are only little differences with respect to 
the best-fit model in the optical spectrum. Although this 
solution cannot be formally excluded a priori, we consider 
it very unlikely, since extinction values so much higher for 
all the intermediate age compared with the youngest pop- 
ulations do not have any physical meaning. 



nous mass. On top of this there is the contribution of a 
strong and heavily extinguished burst of star formation, 
which account for the strong observed emission lines and 
for the intense FIR flux (see Fig.lllJI. but give a negligible 
contribution to the total luminous mass. The extinction 
shows the characteristic age-selective pattern in which the 
youngest populations are more heavily obscured than the 
older ones. In this case the youngest population is com- 
pletely obscured with an E(B-V)~ 5 and cannot be seen in 
the optical spectrum. The old populations with age > 10 9 
yrs have E(B-V)<0.5. 

In order to test the possibility of an alternative star 
formation history, characterized by a much longer star- 
burst duration, we tried a constant star formation rate 
over the last 5 • 10 8 years. The corresponding star forma- 
tion and extinction patterns are shown in Fig. 1131 with 
this constraint, the code attempts to minimize the impact 
of the intermediate age SSPs by increasing their extinc- 
tion. The predicted IR flux is within 3<r from the observed 



7.2. Models for the northern nucleus 

The spectrum of the northern nucleus (spectral regions # 
13 and 14) was analyzed with the same fitting procedure 
described previously. The results are shown in Figs. PHI 
and 1151 the star formation history of this nucleus closely 
resembles that of the southern one, with two main star 
formation episodes and a recent, much weaker, burst. The 
extinction appears to be lower than that of the southern 
nucleus. 

The old population contributes up to 75% in mass and 
the intermediate age population (0.5 Gyr) about 25%. 

Note that there are problems in fitting the observed 
continuum in the range 6000 A to 7000 A, which wc were 
not able to solve, even adopting different extinction laws, 
SSP metallicity, or different SSP combinations. A simi- 
lar problem is found between 3400 and 3800 A, where 
the model systematically overestimates the observed spec- 
trum, but here the observed spectrum has a low S/N ratio. 
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Table 3. Spatial dependence of the star formation rate (M /yr) and extinction (Ay) in Superantennae, as derived 
from the best fit models. Each row refers to a different SSP; ages and durations, in years, are reported in columns 
two and three respectively. The following pairs of columns refers to the numbered galaxy regions and provide the SFR 
and extinction values for each SSP. Values are not shown when the star formation rate is negligible. The bottom three 
rows in each panel give the total stellar mass sampled by the spectrum and the Ay extinction as computed combining 
the three youngest populations (young) and all the 10 SSPs (total). See text for more details. 
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Fig. 14. Best fit to the northern nucleus spectrum. The 
mismatch at the edges of the spectral range analysed is 
found in all other models we have tried. 



7.3. Models for the regions outside the nuclei 

We have obtained best-fit models for the spectra of all 
the individual regions along the line connecting the two 
nuclei. Here both the continuum and line intensities are 
much lower than those from the two main nuclei (Fig.[SJ), 
and the star formation activity correspondingly low (see 
Table |5J. Nevertheless the main characteristic of the star 
formation history of the two nuclei, namely the evidence 
for an old stellar population and for two main episodes of 
SF, one at ~ 10 9 yrs and one ongoing, seems also reflected 
in the spectra of the extranuclear regions. 

8. DISCUSSION 

We discuss here our previous results on the star formation 
history, extinction properties, and the AGN contribution 
in the Superantennae. 

8.1. Spatial-dependent star formation and extinction in 

Superantennae 

Table reports our modelled values of the average star 
formation rate (M Q /yr) along the slit for stellar popu- 
lations at different ages. This quantity has been derived 
from the total stellar mass of the population of a given 
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Fig. 15. SFR and colour excess of the northern nucleus 
best fit model. The star formation history of this nucleus 
closely resembles that of the southern nucleus, with two 
main star formation episodes, forming 75% and 25% of 
the total stellar mass respectively. The intemediate age 
population provides the bulk of the FIR flux. Extinction 
is significately lower than in the southern nucleus. 

age divided by the duration At of the SF episod (third 
column in the table). We have then computed the total 
stellar mass, sampled along the slit, and the correspond- 
ing visual extinction for the various populations. 

Figure ITHI compares the E(B-V) values obtained with 
our global fitting procedure with those derived from the 
Balmer decrement, after correcting H/3 for stellar absorp- 
tion. In the upper panel we show the E(B-V) distribu- 
tion obtained combining the three youngest SSPs, while 
in the lower panel that derived averaging over all ten stel- 
lar populations. In the latter case the Balmer decrement 
technique appears to sistematically overestimate the ex- 
tinction, while in the former there is a better agreement 
at least close to the southern nucleus. This is because of 
the different extinction properties of various populations, 
which the global fitting model can deal with, but not the 
integrated spectrum. Also, by accounting for the contribu- 
tions of the intermediate age stars, the model provides a 
fair estimate of the reddening even in those regions where 
H/3 is in absorption and the Balmer decrement technique 
is inapplicable (e.g. spectra from 9 to 17). 

Table 01 compares the ongoing SFR values obtained 
from our global fitting procedure (c ombining t h e two 
youngest SSPs) with those derived bv iKennicuttl Jl998) 
calibrations: 

SFR{[Ou}) = 1.4 • 1CT 41 x L ([On]) M Q • yr" 1 

SFR (FIR) = 4.5 • 1CT 44 x L (FIR) M Q ■ yr' 1 ( ' 

and eq. El Before computing the SFR, all the measured 
fluxes were corrected for internal absorption using the ex- 



Fig. 16. The spatial dependence of the extinction, as com- 
puted on the basis of the Balmer decrement and through 
the spectrophotometric modelling of the spectra (see table 
13). The Balmer decrement has not been applied in regions 
9 to 17 where H/3 is not detected or appears in absorption. 

tinction values found in Section 3. A value of E(B-V)=0.4, 
based on the results of the spectral synthesis, was assigned 
to those spectra for which H/3 is not detected or is found to 
be in absorption (see Section E}.1.2|) . These relations refer 
to a Salpeter IMF between 0.1 M Q and 100 M Q , while the 
IMF limits in our code are 0.15 and 120 M Q : the factor be- 
tween the two is J Q l to -1 dm/ j Q 15 m~ 1 - 35 dm ~ 1.162. 
The model's SFR of the southern nucleus exceeds by about 
20% the SFR estimate based on the FIR flux, but this is 
within the uncertainties of the calibrations if we considere 
that eq. refers to an average continuous SFR over the 
last 100 Myr, while the best-fit model in Figure El as- 
sumes a very young short burst: the latter needs a higher 
SFR to supply a given FIR luminosity 

The [On]A3727 diagnostics leads to a higher SFR 
than the Ha, with a larger discrepancy in the circum- 
nuclear regions. The [On]A3727 line has however a very 
different spatial profile from [Oiii]A5007, Ha or H/3 (see 
FigEJi indicating a different degree of ionization, possi- 
bly due to a metallicity gradient (see Fritz et al. 2002 
for a detailed discussion). Moreover the SFR obtained 
by the [Qil]A3727 line has a large intrinsic uncertainty 
llCharlot fc LonghettiL l200l|) . A fraction of the star for- 
mation in the southern nucleus and in the circumnuclear 
regions is apparently hidden (SFR from Ha is smaller than 
SFR from FIR even corrected for the AGN contribution) , 
i.e. is not recoverable from the optical spectrum. 

Concerning the northern nucleus, the estimates based 
on the line fluxes and on the spectrophotometric synthe- 
sis model analysis agree fairly well. Apparently this region 
of the galaxy is likely in a post star-burst phase, as sug- 
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gested by the virtual absence of a very young component. 
In this substantial discrepancy between the SFR 

value based on the FIR flux and those based on the other 
methods is apparent. Figure ED indicates that the FIR 
emission of the northern nucleus is produced mainly by a 
5 • 10 s years old population; therefore applying the SFR 
calibration based on FIR leads to a misleading estimate, 
while our modelled SFR value appropriately accounts for 
the ongoing star-formation activity. Because the youngest 
populations are those producing nebular emission lines, 
the Ha and [On] diagnostics are in much better agree- 
ment with our modellistic value. 

Concerning the past star formation history, we find 
evidence for two distinct episodes in this merging sys- 
tem. A recent, still ongoing burst producing the opti- 
cal emission lines and the strong infrared flux, and an 
older one which took place about 1 Gyr ago. Remarkably 
enough, this double-peaked recent star formation history 
appears to be common to all the regions within our aper- 
ture. Ana lysing a sample of int eracting Ultraluminous IR 
Galaxies, iMurphv et al.1 l)200l|) suggest that the merging 
events may evolve through various steps corresponding to 
different encounters, hence ULI RGs may show eviden ce for 
multimodal stellar generations. IColina et a,lJ (|l99lh give 
an estimate of about 1 Gyr for the age of the first en- 
counter, based on the extent and the speed of the gas 
in the tails of the system. Our results are fully consistent 
with this overall picture of the merging process, and relate 
the two main episodes of activity with the first encounter 
and the final merging of the two galaxies. In this context, 
Figure ITol suggest that at the time of the second enounter 
the gas in the northern nucleus was already almost ex- 
hausted. During the first encounter about 30% of the to- 
tal stellar mass was formed, while the remaining 70% had 
already been formed by a pre-existing older generation of 
stars. 

Our model gives very similar ma sses for the two nu- 
clei, consistent with the analysis of IColina et al ] lll99lh 
who, based on dynamical considerations, infer a ratio ~1 
between the masses of the two progenitor galaxies. 

8.2. The contribution of the AGN 

Our analysis of the optical-FIR SED has found evidence 
for a substantial type-2 AGN component in the southern 
nucleus. We have shown, in particular, that the J-L colour 
may be exploited as an estimator of the AGN contribu- 
tion to the MID and FIR emission. For example, while 
the northern nucleus has J-L~2.7 value entirely consistent 
with a galaxy model dominated by evolved stellar popula- 
tions, the southern nucleus shows an excess in the L band, 
J- L~3.9 mag, that unequivocally proves the presence of 
an AGN. Fitting the global SED with a combination of 
an AGN and starburst templates, the contribution of the 
AGN has been constrained to be 40%— 50% of the bolo- 
mctric flux between 8 and 1000 ^m. The uncertainty in 
this estimate is mostly due to uncertainties in the param- 
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*: values obtained assuming E(B — V) = 0.4 
in these regions of the galaxy. 

Table 4. Values of the ongoing star formation rate SFR 
(in M Q /yr) derived from: the main emission lines, after 
correction for extinction and for the AGN contribution, 
and adopting the the Kennicutt's calibration; from the 
far-IR luminosity (Kennicutt's calibration); and from our 
global best-fit model. The latter has been computed com- 
bining the two youngest single stellar populations. All 
SFRs dignostics adopt a Salpeter IMF between 0.1 and 
100 M . 



eter R setting the scale size of the circum-nuclear dust 
distribution, which affects the dust temperature distribu- 
tion. 

The optical spectrum also indicates the presence of 
a Seyfert-2 nucleus. Pernechele et al (2002) provide evi- 
dence for ~2% scattered polarized light around Ha in the 
southern nucleus. Also the analysis of the [Oiii]/H/3 ra- 
tio indicates a non-thermal contribution of ~ 40% at H/3, 
consistent with that required by the MID IR excess. 

In X-rays, Pappa et al. (2000) showed that 
Superantennae is characterized by a very hard 2—10 keV 
spectrum with Lx — 2x 10 42 erg/s, and Braito et al. (2002) 
found evidence for a high-equivalent width Fel Ka line 
typical of heavily obscured type-2 AGNs. 

At radio wavelengths, the contribution of the AGN 
may be evaluated by comparing the observed flux with 
that expected from a starburst galaxy on the basis of the 
FIR/radio correlation (e.g. Sanders & Mirabel, 1996): 



q = log 



F FIR /(3.75 x 10 12 Hz) 
S„(1.49GHz)/(W m- 2 Hz _1 ) 



2.35 ±0.15, 



(5) 



where F F1R = 1.26 10- 14 (2.58 S 60Mm [Jy] + S 100im [Jy]) 
W m~ 2 . Wright et al.(1994) report a radio emission of 
S„ ~ 102mJy at 4.85 GHz. Assuming a radio slope 
a#=0.8, typical of starforming galaxies, the ratio in 
eq. © becomes q~2.71. We find instead from Table [21 
^4 85 GHz —1-83, significantly lower than the expected 
value for a normal starforming galaxy. Conversely, assum- 
ing that the starburst is responsible of ~60% of the total 
infrared flux from Superantennae, and from eq. JHJ), we 
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expect for a pure starburst component a 4.85 GHz flux of 
7.8 mjy, more than one order of magnitude fainter than 
observed. So the AGN dominates the radio emission. 

Finally, the two merging/interacting systems appear 
to have almost the same stellar mass (~6xl0 10 M Q 
and ~3.4x 10 10 Mq). From the correlation of the central 
Massive Dark Object (MDO) mass in galaxy nuclei with 
the bulge mass (Magorrian et al. 1998), we would then 
expect that the two nuclei host similar-mass MDO's. The 
fact that only the southern nucleus is site of AGN ac- 
tivity is consistent with our discussed evidence that the 
northern nucleus has underwent a past activity episode 
(possibly involving both a starburst and an AGN), whose 
feedback and/or gas exhaustion have possibly quenched 
out any subsequent activity. 

9. CONCLUSIONS 

We have presented a spatially-resolved study of the star 
formation and extinction properties of the Ultra Luminous 
InfraRed Galaxy IRAS 19254-7245 (The Superantennae) 
along the line connecting the two nuclei. 

We have combined new optical spectroscopic obser- 
vations with detailed broad band spectral energy distri- 
butions built from published photometry. We have first 
attempted to quantify the presence of an obscured AGN 
- for which there are evidences from both X-ray, opti- 
cal and radio observations - and its contribution to the 
source spectrum: we have found that the J-L colour may 
be a powerful diagnostic tool. By means of an elaborate 
fit of the whole SEDs, from the optical to the FIR, we 
have inferred that about 40-50% of the luminosity of the 
southern nucleus in the infrared is due to an obscured 
AGN. The northern nucleus instead does not show any 
evidence of an AGN. 

After having accounted for the AGN contribution, the 
combination of the optical spectral data with infrared ob- 
servations allows to disentangle between extinction, age 
and star-formation effects: our spatially- and temporally- 
resolved analysis is able to largely resolve all the degen- 
eracies otherwise present in the integrated spectra. The 
data on Superantennae are consistent with a recent bi- 
modal star formation history. The first episode of stellar 
formation possibly corresponds to a first collision between 
the two progenitor galaxies, about 1 Gyr ago, when a sig- 
nificant fraction of the gas mass was converted in stars 
(about 30% of the present total mass of the system) . The 
merging onto the current configuration is due to a more 
recent episode, adding a negligible fraction (few %) of the 
total stellar mass. 

We show that the northern nucleus has been character- 
ized in the past by a slightly older and less intense burst 
of star formation. 

The radio emission by Superantennae, which is more 
than one order of magnitude larger than expected from 
the starburst if we use the FIR/radio correlation, con- 
firms the presence of an AGN. Having said that, the radio 
flux cannot be used to estimate the AGN contribution 



to the bolometric emission. Our analysis also emphasizes 
that the technique used to estimate the redshift for dusty 
starbursts and based on the radio to mm flux ratio may 
provide highly unreliable results if even a minor AGN com- 
ponent is present in the system. 
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